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1, ?i:r..>rgj ! [tacwsg ton ? 

The problem Is to track Che motion of part of a field of view. Let, 
os assume that the scene Is a twa-diaensicnal picture in a plana perpen¬ 
dicular to the roll axis, (These Sir.plifyi.ng AS sumptions t of course, stc 
a main problem is estimating how thE system works in real life). So we 
tan think of the picture as a function f(x,y) in seme plane- 

How suppose that at time t Q thE scene Is f 0 (x,y) and at some time 
inter it has icved, and is f^(sL^y). Suppose also that the seen® has not 
changed, hut has only teen moved rigidly in the plane. Then an elegant 
ma th enat i c al. way to as ti mat a this (notion Is to compute (he cross-coir Elation 
of the original and current picture. First let us review a basic simple 
mathematical fact. Given any function £(*) nod OOv displacement A, it is 
true that 

UP 3D 

/ f(K> f(n] £ ; f U) f (x +4) 

— tc- _CD 

(when the expressions are meaningful) and we have. > rather than = 
except under the most peculiar conditions, 1►e,, when the pattern is 
perfectly periodic, (This fact is a consequence of the ttiAngla inequality 
a^ + h^ > c^, slightly generalised). In fact, if one considers 

(i) - /*(*) f<K V A) 

■one always gets a graph likfl 



with its maximum at 0 and symmetrical 00, both sices, How, consider soffe 
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small number d t called the "delay 11 ', and consider the formula 

fi(h) = n f(x + d) - ft* - d) ] ftK + h) 

If we look at the values of this for various Values of h, we find: 

g{d) = ff(x +■ d) ffr +■ d> - /f (k - d> £(k + d) = +(o> - tf(2d) *0 

e(o) - Itu + d) 100 - m« - d> f00 = ^fd> - $(d > - o 

e(-d> " / f(x + d) f(x - d] - /f (x ~ d> f(x - d) = $(Zd) - tf.(o) <0 

<» W 

{because f£ (x ■+ a) f £x +■ a + b) = /f (x) f (x + b) in general) . 

- -40 

and assuming everything is nice and smooth, we nan expect CO get a Curve 

like this: 



yielding a "discrimination curvet" which is just what we want, because in 
the ’"linear" range it computes the displacement of f(x) from f (xH-h) for us 5 
How it ia of practical importance for us to know how large is the "linear” 
range hecause 

(1) it tells us how large a step a device could cope with and 

(2j it incidentally tells ua how accurate or coarse the computation 
of the integrals must be - this bears on the size of optical 
slit - or whatever is used in practice. 

Hew we can find a let about this by (theoretically) considering a very 
snail displacement d - in fact, if d is allowed to go to aero* we get a 
sort of differential analysis? we divide by d to keep things from blowing 
up h and we ge t : 
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llm 

d+o 


/[£(* + d} - f (x - d)] f (jt + h) 


“ a 


f(* + d) - ftM - d> 


f (x + h) 


2d 


= 2 / f 1 (ac) ft* + Ei} 

- 2 +' (h) 


by soot analysis involving de(h) 

dh 


and Che fact that $(h) ■* 


f f' (x)f(x +■ h} + f r (x)f(x + b)j. 


This a ays that as the delay d goes to zero, the n dia crimination curve 1 ' 
approaches the derivative o£ the auto-correlation function *(h) 



*n<3 this suggests that there is no serious risk if we take i to be of the 
order of % the width of the Drain peak of the a u-to-cor relation £l netion. 

In terms of tb.E picture-function £(x,y), this width is a measure of the 
"grain" of the scene-texture; it is probably the most important ("first- 
order") parameter for characterizing scenes in terms suitable for 
correlation-based ins Lrucents , 

What is the intuitive significance of the tracking function? 

jTf'fa) £(k + h) = £ r (h) 

Two simple examples show who E is happening i 
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1. =_i. 

r 00 - 


so if h is positive* or negative we get 




and here, for a square wave t we .get simply the sian of the dlapLacsiretit l 
on the other hand h 



b positive: 


h negative. 



and we get something mere proportional to the artOvnt of the displacement. 
Under the assumption that the lengths of regions of positive and negative 
derivatives ate generally e^al ( we can approximate this "analogue 17 
situation by differentiating and ’'clipping 11 the results] the error in 
this depends very much on the situation and sometimes the "digitised" 
result will be better, usually worse, titan the analogue. 

£. P olar coordinates and astial ttacking . 

Now 1st us consider a more practical, situation; w£ M ns id ct the 
prebles of tracking, still with the same assumptions, but assuming the 
scene i 3 on a sphere sc that we want to compute pitch, yaw* and roll.. 

Now let us replace L(jt s yO by Ffr*0) and consider a polar integral 



"■ 2rr 

s E(Tf) ; Q(r,B) 4(B) d &dr 

* o 

wbor* X. and J ace 'hardware" fu.c.c,tions ue have to consider. Suppose 
first that J C-S} = I and K(if) falls cff fast enough to make the Integral 
converge (that £ 3 , a "finite radial angle shaded slit." the slit width 
is concealed in Che real of the mathometits) T Then If the scene is 
simply rotated, wv use for Qir,^) 

qtr, 9) - F L (r, 6 ) F(r,0 +k) 

and this gives ug an egtimate of the phase *££ , Just as in tEne first 
section► 

If the scene suffers a translation , in pitch or yaw* the effect 
on this averages Out to aero because we get opposite effects in each 
hecisphere: 



but this is a statistical cancellation and is an error that must be 
reckoned out Thus there will he "cress talk" if the picture is not 


symmetric „ in a suitable average sense. 

Now* to get the yaw-Component, we use J(£> ■ sin5 and get the effects 
above to add t 



+■ (advance) 
- (retard) 


sin = + 
4 



In this case the contribution of the pitch component averages to aero 
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he.cau.ae. for it. we get cos (S) ■ J(0); 

+ [advance) ccs(^) = ■+ 

7T 

- [retard) cos * - 

acid these cancel. 

Thus it is possible to resolve the three Components roll (J(S) 31 1) } 
pitch, (J(0) = cos $) , and yaw (J (G) E sin 0) by changing the 1 Vernel 1 
of the integral. One Could even approximate this digitally by using 


i 



A3 approwlitiations, hut this would increase the "cross talk*" 

3-. the As Son? t Ions 

It remains to discuss the assumptions; 

1 H Stereo changes mean that throughout, we never really have 
f(x -&■ h) 

but w£ always have a new picture 
g(n + h) f(x + h) 4- error 

2 , It i-s safe to say chat because the situation is real, there is 
no trouble with die analytic mathematics, One has to worry only 
about the photon and othar fihisa involved in empirical evaluation 
of the inCagrals and derivatives, 

3. Although che auto-cotrelation function is always symmetrical 

/f(x> fU + h) « /ffci) £[x - h) 

this Is only approximately true for a nev scene 
/g[jO f(x +■ h) ^ /&(?;) £{k - h> 
anu a systematic drift during a changing scene could CAUSA A 


cumulative error. 
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4* .The most serious problem is: 7a what extent will a theoretical 
calculation based on a good smooth* convex maximum of auto- 
correlation function be realistic? It tflLL b* Valuable in 
calculating, the underlying, basic random^variable drift of the 
System.* but this theory has to be supplemented by a special- 
case analysis of the worst "deterministic' patterns known or 
suspected to cause trouble. 

We should do some experiments using correlation on very simple 
arrays. Evan, say * Lfi points around a circle should yield interesting 
results - 



